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1 Introduction and motivation 



One of the most interesting experimental results in cosmology in the last years has been the 
measure of the acceleration of the universe from the observation of high-redshift supernovae 
0]. The combination of this result with the measure of the cosmic microwave background 
(CMB) anisotropies favors a universe whose energy density is very close to the critical 
one, fitotai — 1) and with a large dark energy component, fidark — 0.7. This realization that 
the energy content of the universe might be dominated by some kind of dark energy has 
revived the interest in cosmological models where the late time expansion of the universe is 
governed by some kind of vacuum energy, either in the form of a nonvanishing cosmological 
constant or a quintessence field (see 0, |J] for comprehensive reviews). 

Recently, braneworld scenarios have emerged as a very appealing alternative to standard 
dimensional reduction in higher-dimensional theories. The general idea is the following. The 
observable universe is a submanifold embedded in a higher-dimensional space-time. By some 
physical mechanism matter fields are confined to this submanifold. Although gravity can 
propagate into the bulk, the Kaluza-Klein excitations of the four-dimensional graviton can 
be decoupled at low energies either by considering finite but "large" extra dimensions f5[ or 
a negative cosmological constant in the bulk |J. Therefore, at energies below the mass gap 
of these Kaluza-Klein modes, physics on the brane is well described by the matter degrees 
of freedom confined to the brane and standard gravity in four dimensions. One of the most 
interesting features of these models is that the energy scale at which the five- dimensional 
effects are important can be relatively low, even of the order of the TeV, and therefore in 
principle accessible to future high-energy experiments. 

It has been realized that the braneworld paradigm in its cosmological application leads 
to corrections to the four- dimensional Einstein equations on the brane |7], [|, ||. These 
departures from General Relativity, and their imprints in braneworld cosmological models, 
might be used to test the scenario against cosmological observations. In the context of 
Friedmann-Robertson- Walker cosmologies, the braneworld hypothesis results in the presence 
of two extra terms on the right-hand side of the four-dimensional Friedmann equation: a 
dark-radiation contribution and a quadratic term in the energy density that couples directly 
to the five-dimensional Planck scale. While the first one might have any effect at the present 
stage of the history of the universe, the second one is expected to be important only in the 
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very early universe where the energy density was very large. 

Quintessence models have also been studied in the context of brane cosmologies [10, |TT|| . 
With a few exceptions |12|, in most of the models considered a four-dimensional decaying 
cosmological constant is implemented by introducing an additional scalar field in four di- 
mensions. By a suitable choice of potential its contribution to the energy density becomes 
dominant only at present. 

In this paper, we study the possibility within a braneworld scenario that an effective 
four- dimensional quintessence is induced by a dilaton field propagating in the bulk. In our 
model the cosmological evolution in four dimensions reflects the movement of the braneworld 
in a static dilatonic bulk space-time in which Z2-reflection symmetry around the location 
of the brane can be broken. We find that the bulk dilaton induces a term in the effective 
four- dimensional Einstein equations that acts as a decaying cosmological constant. Hence, 
in this sense quintessence occurs naturally in this class of models. For universes filled with 
matter and radiation, and for a certain range of the parameters defining the model, this 
vacuum energy actually tracks the matter or radiation energy density and dominates at 
late times, driving an accelerating expansion of the universe. Using constraints from big- 
bang nucleosynthesis and the present value of Newton's constant allows us to constrain the 
parameters of our class of models. In particular, we find a lower bound of 10 5 GeV for the 
five-dimensional Planck mass. 

In the solutions showing quintessencelike behavior the universe first undergoes decelerated 
expansion and then at late times enters into a regime of accelerated expansion due to the 
effects of the bulk dilaton. In this case, the observed isotropy of the universe has to be 
assumed. However, as will be seen, the class of models under study also allows for inflationary 
solutions at early times. In this type of solution the universe starts in a state of accelerated 
expansion and then later undergoes decelerated expansion. The problem in this kind of 
models, however, is not only the lack of acceleration in the expansion at late times, but the 
fact that in this asymptotic regime the decelerated expansion is dominated by the p 2 term 
in the Friedmann equation. 

Much of the recent interest in brane cosmologies was triggered by Hofava-Witten type 
models |13j where the brane is confined to lie on the fixed points of an S 1 /Z 2 orbifold. Thus 
usually a Z2 symmetry in the bulk direction is assumed. However here we will consider more 
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general situations in which the bulk space-time is not symmetric under reflection around the 
position of the brane. In our context, this situation describes braneworlds that emerge as 
domain walls separating two regions of space-time with different values of the cosmological 
constant. As we will see, this asymmetry leaves its imprints in the cosmological dynamics on 
the brane at late times, for example through corrections to the four-dimensional Newton's 
constant depending on the difference between the cosmological constants on the two sides of 
the brane. 

The plan of the paper is as follows. In Section 2 we study the dynamics of moving 
branewolds in the background of five- dimensional Einstein-Liouville solutions. Section 3 will 
be devoted to the study of the cosmological dynamics of moving brane universes with 
reflection symmetry, whereas in Section 4 we will analyze brane models where this reflection 
symmetry is broken. In Section 5 we will proceed to study the phenomenology of these models 
by using both current observational data and bounds derived from big-bang nucleosynthesis. 
Finally, in Section 6, we will summarize our results. 

2 Moving branes in a static bulk geometry 

In Ref. [|14| a broad class of braneworld cosmologies coupled to a dilaton field in the bulk 



was studied by using a solution generating algorithm that enables the construction of five- 



dimensional Einstein-Liouville cosmologies (cf. Assuming Z 2 -reflection symmetry 

around the fixed location of the brane, it was found that the solutions studied enjoy a new 
version of the self-tuning mechanism of in which the total cosmological constant in 



four- dimensions is canceled independently of the value of the bulk cosmological constant 
and the energy density on the brane. It was also found there that inflation in Friedmann- 
Robert son- Walker and Bianchi-I models belonging to the family of solutions studied cannot 
be achieved for any value of the amplitude and slope of the bulk Liouville potential. 

Here we are going to use these bulk solutions in a different physical setup. Our aim is 
the construction of realistic four- dimensional brane cosmologies in which the cosmological 
expansion accelerates at late times. For this, even if we restrict the bulk solutions to those 



studied in D14J , we will have to relax some of the conditions imposed on the braneworlds 
constructed there. In particular, here we will not consider the braneworld to have a fixed 
location on the bulk, but we will take into account the degree of freedom associated with 
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its collective motion in the five-dimensional space-time. In addition, we will also consider 



situations in which Z 2 -reflection symmetry is broken (cf. flT5| , jig , As we will see in the 
following, the features of the resulting brane cosmologies are quite different from the ones 



studied in [14 



2.1 Description of the model 

We will consider solutions to the five-dimensional bulk action 
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(2.1) 



Here we are assuming that the five-dimensional Planck scale k 5 determines also the scale of 
the five-dimensional scalar field, which is dimensionless here. The only other independent 
energy scales entering into the action (|2.1| ) are provided by A, which has dimensions of 



(mass) 2 , and the energy scale of the four-dimensional dynamics. In the following, we will 
consider that the matter on the brane is described by a perfect fluid with barotropic equation 
of state p = wp (with \w\ < 1), so the matter Lagrangian is written in terms of the energy- 
density as ^matter = —p- In |T4]] , solutions to the equations of motion in the bulk derived 
from p.lj ) were constructed in terms of the functions Q(x) and h^ u (x) in which the metric 
takes the form[J 
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(2.2) 



with the scalar field given by 



4>(x, r) 
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Q(x) + fclogr. 



(2.3) 



The metric function h^ u is a solution to the four-dimensional Einstein equations coupled to 
a scalar field, R^ v = d^d^ with ip(x) = y6Q(x). The amplitude of the Liouville potential 
function of k and £ has the following form 

1 



A = -(k 2 -u)e. 



(2.4) 



4 These solutions are related to the ones studied in [Q by a fc-dependent coordinate redefinition r 2 
exp (a 2 x) an d a shift in the dilaton field. 
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The parameter £ has dimensions of mass and determines the scale of the cosmological con- 
stant. 

Here we will study the motion of a braneworld in a bulk geometry described by the line 
element ( |2.2|) , possibly with different values of £ on the two sides of the brane. In order to 
simplify our analysis, we consider Q(x) =0 in Eq. ( |2.3|) and the bulk space-time to be static. 
This implies in particular that is a static Ricci-flat four-dimensional metric. Hence the 
metric takes the form 

ds 2 5D = l r K fc2 - 3 )dr 2 + r 2 h^dx^dx", (2.5) 

while the dilaton is 

4> = k log r. (2.6) 

In addition, the four- dimensional metric is given by 

h tiv dx il dx v = -dt 2 + hij(x)dx l dx j . (2.7) 

Let us then consider a single braneworld embedded in the five- dimensional space-time 
described by Eq. ( |2.5| ) and moving along the trajectory 

t — t(r), r = R(r), x, y, z = constant, (2.8) 

where r is the proper time. For a five- dimensional static metric of the form ds 2 = Q^dr 2 + 
Qoodt 2 + Qijdx l dx\ the nonvanishing components of the 5-velocity V A and the vector normal 
to the brane n A can be computed using the relations VaV a = — 1, Vxn A = 0, with the result 



d / ^55 5 / 1 + G^R 2 /„ n x 

V — WOO V ^55 

Here and in the rest of the paper the dots indicate differentiation with respect to proper 
time t. 

On the other hand, taking an orthonormal frame {Vm > n M, e ^ii}> ^ ne components of the 
extrinsic curvature can be easily computed using the existence of the timelike Killing vector 
of the static bulk space-time (cf. 0, f21|j): 

1 d 



Kl 



\T I G^Jl + Q^R 2 



(2.10) 
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k) = ^n M g ik d M g k3 . (2.11) 
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Note that K TT is given by the orthogonal component of the proper acceleration of the brane. 

Looking at Eqs. ( |2.5|) and (|2.7|) , we see that, when expressed in proper time, the induced 
metric on the brane takes the form 

dsl B = -dr 2 + R 2 (T)h ij (x)dx i dx j . (2.12) 

Since the metric (|2.7|) has to be a static vacuum solution to Einstein equations, in the 
following we will focus on the simplest possibility where hij = 5ij. In this case the line 
element ( |2.12| ) describes a flat Friedmann-Robertson- Walker universe. 



2.2 Junction conditions and dynamical equations 

As explained above, we will consider that the constant £ determining the amplitude of the 
Liouville potential takes in general different values ^± on the two sides of the brane (a related 
proposal has been studied in |16|)- This means that the bulk geometry is described by the 
line elements 



d 4D 



d 4D 



S + 

= 2- r f (* 2 - 3 ) dr 2 + r 2 [-N(t) 2 dt 2 + Sijdx'dx'], r < R(r). (2,3, 

Following |16] we have introduced a shift function N(t) in order to ensure the continuity of 
the metric across the brane. Notice that the bulk field equations are satisfied since, given the 
static character of the metric, the shift function can be locally reabsorbed in a redefinition 
of time. A simple calculation shows that the metric is continuous across r = R{t) by taking 

m = (y e + a^m (2 ' 14) 

On the other hand, the dilaton is continuous across the location of the brane. 

The Israel junction conditions, together with the matching for the dilaton across the 
brane, will determine the dynamics of the scale factor R(t). As for the metric junction 
equations we find 

[K, u ] = K+ - K- u = ^ [-A(0) + g^T aP ] g, v - «1 V , (2.15) 

where g^ v is the induced metric, t^ v is the energy- momentum tensor on the brane and A(0) 
is the brane tension. Indices are raised and lowered using g^ u . Since we are assuming 
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that the matter on the brane is a perfect fluid the energy-momentum tensor is given by 
= diag(— p, p, p, p). Using the barotropic equation of state p = wp we find from the 
matching condition for K^, 
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(2.16) 



while from the jump of at r = R(t) we get 



3 A; 2 
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(2.17) 



In addition to this, we have also to impose the matching condition for the dilaton derivative 

(2.18) 



across the brane 21 



[n A d A (f\ =k 2 5 [A'(0)-&<7>"V] 



which, upon substitution, reads 
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K 2 5 [\'(<P)-(3w-l)bp]. (2.19) 



The particular case in which the bulk space-time has Z2-reflection symmetry across the 
position of the brane is recovered by setting £ + = = £ and taking the outward normals 
on the two sides of the brane with opposite signs. In Eqs. ( |2.16| )-( P.19| ), this amounts to 
replacing the minus sign between the two terms on the left hand side of these equations by 
a plus sign. 

Equations (|2.16| )- (|2.19| ) derived from the junction conditions determine the dynamics of 
the four-dimensional scale factor R(t) in terms of the brane tension \(4>) and the energy 
density p. Consistency of Eqs. ( |2.17| ) and ( |2.19| ) implies the following differential equation 
relating the brane tension with the energy density: 



A'(0) + |a(0 



(3w-l)b 



k 



^2.20) 
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It is remarkable that, because of the junction condition for the dilaton, the matter energy 
density on the brane actually mixes with the brane tension. 

It is important to stress that any solution of Eq. (|2.17| ) will automatically satisfy the 



second order differential equation derived from Eq. fl2.16|) , provided the brane tension X(4>) 
and the density p satisfy fl2.20|) and ( |2.22| ). Thus, as it is the case in standard cosmology, the 



dynamics of the brane universe is solely determined by the first-order differential equation 

The analysis of the case k = requires some care. Since for vanishing k the dilaton field 
exactly vanishes for the family of solutions studied, the dilaton matching equation (|2.19| ) 



is not imposed and the only dynamical equation is Eq. (|2.17|) with constant brane tension 
X(4>) = Xq. In this case, therefore, there is no mixing between the vacuum energy and the 
matter on the brane. 

2.3 Energy balance on the brane 

In order to get the form of the brane tension A(0) from Eq. ( |2.20|) it would be necessary 
to know how the energy density depends on the scale factor (or equivalently, for k ^ 0, on 
the dilaton field). The conservation equations for the energy-momentum on the brane can 
be derived from the Codazzi equation. Actually, as explained in M, in the presence of a 
dilaton field there is an interchange of energy-momentum between the brane and the bulk. 
This implies the existence of a nonvanishing right-hand side in the conservation equation 

9 + 3|(p + V) = -(3p - p) b R^. (2.21) 

We see that the source term in the continuity equation results from the conformal coupling 
between the matter-energy density on the brane and the bulk dilaton, measured by b. 

Using the equation of state of the fluid and the expression of the dilaton field in terms 
of the scale factor, Eq. (|2.6Q , we find the following differential equation for p(R): 



R% + mp = 0, with m = 3(w + l) + (3w-l)kb. (2.22) 
dR 

This implies that the energy density depends on the scale factor and the dilaton field as 

p(R)=p R~ m =s- p(0) = p oe -^, (2.23) 
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where the second equation is only valid for k 0. By comparing these expressions with 
the ones encountered in ordinary cosmology, we find that the usual evolution of the energy 
density with the scale factor is corrected by a term that depends on the conformal coupling 
b between the dilaton and the matter fields. Equation (|2.23|) can be physically interpreted 
by considering that, for an observer on the brane the source term in Eq. (|2.25|) leads to an 
effective barotropic index w e ff governing the evolution of the energy density with the scale 
factor p(R) ~ R-^«+ l ) , where 

w efi = w + ( w - \ ) kb. (2.24) 



3, 

In Section 5, we will see that for phenomenologically viable models k 2 < 10~ 2 . Since it 
is natural to expect the coupling b to be of the same order as k, the difference between w 
and w e ff would be also of order 1CT 2 . There are, however, two particular instances in which 
the correction to w in Eq. ( |2.24j ) vanishes for any value of b. The first one is when k = 0, 
which corresponds to switching the dilaton off, so the right hand side of Eq. (|2.21|) is zero 
and the usual continuity equation for p(R) is recovered. The second one corresponds to the 
case in which we consider a radiation fluid on the brane, w = |, where again the right-hand 
side of Eq. ( |2.21|) vanishes. This is a consequence of the classical conformal invariance of 
radiation, which is therefore insensitive to the conformal coupling to the dilaton. As a result, 
when the matter on the brane is dominated by a gas of massless particles, the temperature 
evolves with the scale factor in the same way as in radiation-dominated standard cosmology, 
T(R) ~ R- 1 . 

By plugging the expression for p(0) into the right hand side of Eq. fl2.20|) we find a 
differential equation for the brane tension. When k 2 ^ 0, 3m, the solution is given by 

9(w + 1)' 



A(0) = A e 3 



1 + 



k 2 — 3m 



poe"^. (2.25) 



where Ao is an integration constant. In the particular case when k 2 = 3m the solution takes 
the form 



A(0) 



k 2 =3m 



3 

A - t(w + !)Po0 



k 



e"3 . (2.26) 



Finally, if k = the brane tension is constant and does not mix with the matter-energy 
density, 

A(0) fe=o = A . (2.27) 
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In the following, unless stated otherwise, we will only study the case when k 2 ^ 3m. 

As pointed out above, it is interesting that the expressions (|2~25l) and (ggg) for the brane 
tension depend on the matter energy density through the constant po- This mixing between 
matter and vacuum energies on the brane can be seen to be dilaton-mediated, in the sense 
that it is produced by the junction conditions of the dilaton on the brane, which couples to 
both the brane vacuum energy (through its dependence on 0) and the matter fields in the 



action ( |2.1| ). Only when k = is the dilaton switched off and there is no mixing between 
the matter and the vacuum energy on the brane. 



Moving braneworlds in a Z2- symmetric bulk space- 
time 



In this section, we are going to analyze the cosmological dynamics of a four- dimensional 
braneworld moving on a bulk space-time with Z 2 -reflection symmetry around the location of 
the brane. This case is of some interest, since it represents a deformation of the cosmological 
models studied in [TJJ] by taking into account the degrees of freedom associated with the 
movement of the brane in the bulk. 

The dynamical equations for the scale factor are obtained from Eq. Q2.171 ) by taking 
£ + = = £ and reversing the minus sign on the left hand side. This results in 

= -£ 2 iH fc2 + |(A + p) 2 . (3.1) 
35 

As explained in the previous section, there is a mixing between the brane tension A(0) and 
the energy density due to the mediating dilaton. This implies that the expression of the 
tension as a function of the dilaton field ( |2.25| ) contains a term depending on p. If we collect 
all the p-dependent terms in Eq. ( |3.1| ) together we can write the Friedmann-like equation 
(k 2 ^ 3m) 





36^° 



36£ 2 )ir§ fc2 + ^(±A iT 



9(10 + 1; 



k 2 — 3m 



p(R) 



+ 



^5 
36 



9(w + I) 
k 2 — 3m 



1 2 



p{R)\ 



(3.2) 



where the ± signs correspond to k 2 < 3m and k 2 > 3m, respectively. Let us also remark here 
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that in order for the (time-dependent) Newton constant to be positive, we have to restrict 
A > when k 2 < 3m and A < if k 2 > 3m. 

When w > — 1 we can transform Eq. ( p.2|) into the usual Friedmann equation in brane 
cosmology (cf. p2], ^3|) by defining the effective energy density in four dimensions through 



p(R), (3.3) 



the rescaling, 



k 2 — 3m 

which amounts to a redefinition of the (positive) integration constant po in Eq. ( [2.23 ). In 



physical terms, this means that the energy density measured by a four- dimensional observer 
by looking at the rate of expansion differs from the "bare" value that appears in the action 
( [2.2|) by a positive multiplicative constant. This effective energy density satisfies the same 
continuity equation as p{R) and scales in the same way with the scale factor. We will identify 
it with the energy density driving the cosmological evolution of our brane universe. 

In order to simplify the notation in the following, we will drop the subscript "eff" and 
denote the effective four-dimensional energy density by p(R). Doing so we finally arrive at 
the Friedmann equation on the brane, 

2 4 

= yG w (i?)p Q ( J R) + ^-G N (R)p(R) + |p(i?) 2 . (3.4) 

We see that, in addition to the usual terms linear and quadratic in p(R), there is also a 
"quintessence" energy density given by 

M«) = ^IT^V ^- (3-5) 




As we did in Eq. ( |2.24|) , we can define the effective barotropic index associated with the 



quintessence term that determines the scaling of pq with the scale factor as 

w Q = ^-l. (3.6) 

Thus, for k 2 < 6, it satisfies — 1 < wq < — | and represents a real quintessence contribution 
to the Friedmann equation. In addition to this, the four- dimensional effective gravitational 
coupling in our braneworlds depends on the scale factor according to 

G N (R) = ^irl (3.7) 
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In the case w — — 1 and k 2 7^ 0, the Friedmann equation is retrieved directly from Eq. 
72T) , with the result 

. \ 2 

= - 36 ^)^" §fc2 = Y^(J2)pg(i2), (3.8) 

so the dependence on the matter energy density disappears and the dynamics of the universe 
is controlled only by the induced "quintessence" energy density. When w — — 1 and k — 0, 
there is an explicit dependence on po and the corresponding Friedmann equation is given by 
( p.8|) with the replacement A — > A + Po- 

In order to study the solutions of Eq. ( |3.4j ), it is convenient to rewrite it as the equation 
of motion of a particle moving with zero total energy in a potential, 

B 2 + V(R) = 0, (3.9) 

where 

V(R) = e 2 i?-i (A;2 - 3) - ^ \\\ \R-^ k2 ~ 3) + p R~ {m - l) ] 2 ■ (3.10) 

36 L 

Inflationary solutions on the brane are those with R > 0, R > 0. Using the equations of 
motion (|3.9T), we find that the conditions for inflation reduce to 



V'(R) < 0, (3.11) 

where the prime denotes differentiation with respect to R. 

By looking at Eq. ( |3.10| ), we see that if m > 1, the brane energy density tends to stop 
inflation, whereas the bulk and brane vacuum energies compete with each other. When 
k 2 < 3, the brane vacuum energy contributes to accelerate the expansion, whereas the 
bulk cosmo logical constant tends to decelerate it. If k 2 > 3 the situation is reversed. In 
the following we will study more explicitly the main features of the cosmological solutions 
obtained in the cases when the matter on the brane is dominated by either vacuum energy 
(w = — 1) or a gas of ultrarelativistic particles (w = |). 

3.1 Vacuum branes (w = —1) 

Let us first consider the case in which the matter on the brane contributes to the energy- 
momentum tensor as a vacuum energy with the equation of state p = —p. When k 7^ 0, we 
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find from Eq. ( |3.8| ) that all dependence on p disappears from the Friedmann equation and 
the potential V(R) simplifies to 

V{R)=(?-£)tyirW-*>. (3.12) 

In order to have solutions to Eq. ( |3.9| ) we must demand k^X^ > 36£ 2 . When the inequality 
is saturated (fi^A 2 , = 36£ 2 ) the brane remains at a fixed position R(t) = constant and we 
retrieve Minkowski space-time on the brane for any value of k. The brane tension is given 
by 

A(0) = ±§e-^- Po e 4 ^ (3.13) 
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and the energy-momentum tensor on the brane is = —poe^r)^ (cf. [fL4| ). Actually, these 
are the self-tuning solutions of Ref. [ 19|| which preserve four-dimensional Poincare invariance 



for any value of the cosmological constant in the bulk. 

If, on the other hand, fi^A 2 , > 36£ 2 we find that V'(R) < for k 2 < 3 and hence inflation 
takes place. By solving explicitly the equations of motion ( |3.9| ) with the potential (|3.12|) , we 
find that when k 2 ^ 0, the scale factor scales with the proper time as 



3 



i?(r)~r^. (3.14) 

Thus, for < k 2 < 3, the four- dimensional universe undergoes power-law inflation. 

On the other hand, if k = 0, the brane tension is constant A(0) = Ao and solutions to 
Eq. (|3.9|) will exist provided ftf(Ao + po) 2 > 36£ 2 . When this inequality is saturated, we find 
the static Randall-Sundrum model ||. Otherwise, the four-dimensional metric is not static 
and we have de Sitter inflation, 



R(r) ~ exp | ^Aeff ) , (3.15) 

i.e., in this case the expansion is driven by an effective four-dimensional cosmological constant 
Aefr = f!(Ao + po) 2 -3£ 2 >0. 

3.2 Radiation branes (w = |) 

The second physically interesting case is the one in which the matter on the branes is a 
perfect fluid with equation of state p = |p. The first interesting thing to notice in this case 
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is that m = 4 and all the equations become independent of the parameter b measuring the 
coupling between the dilaton and the matter on the brane. This is a consequence of the 
classical conformal invariance of a radiation perfect fluid. Now the potential V(R) is given 



with A > for k 2 < 12 and A < if k 2 > 12. 

In order to extract the qualitative behavior of the scale factor, we can study the asymp- 
totic limits for large and small scale factors. When k 2 < 12, we find that 



"Large" universes will result only if we demand k^Aq > 36£ 2 . In this case, we observe that 
for large values of R the evolution of the scale factor is identical to the one found for vacuum 
branes. Thus, for k = 0, the universe undergoes de Sitter inflation while for < k 2 < 3 
inflation occurs with the scale factor scaling as a power of the proper time. When k 2 = 3, 
we find coasting behavior of the scale factor and for 3 < k 2 < 12 the expansion of the 
four- dimensional universe is decelerated. 

On the other hand, when R — ► the evolution is always decelerated with R(t) ~ for 
all < k 2 < 12. This is the expected evolution of the scale factor of a radiation universe 
whose dynamics is dominated by the p 2 term in Eq. ( |3.4|) . It is interesting to note that 
the evolution of the universe in this regime is completely dominated by the matter energy 
density, while for large R the expansion rate is determined by a competition between the 
bulk cosmological constant (encoded in £) and the brane vacuum energy as measured by Ao- 

Altogether, we find that for a radiation filled braneworld in a bulk dilaton potential with 
k 2 < 12, the cosmological evolution in four dimensions begins with a p 2 -dominated phase. 
After a period of decelerated expansion resulting in a decrease of the density, the universe will 
enter a regime in which expansion is driven by the term proportional to p in the Friedmann 
equation. This "standard" expansion phase is terminated at some time in which the matter 
energy density decreases below the contribution of the energy density from the quintessence 
term. At this moment, the universe will start its accelerated expansion. To illustrate the 
point, in Fig. 1 we have plotted the potential V(R) and the scale factor for a particular 



by 



V(R) =£ 2 iH (fc2 " 3) 




\ \R-^ 2 -V + p R-* 



(3.16) 




(3.17) 
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model in this class of solutions. In Section 5 we will study this type of brane universe as 
quintessence models for the observed acceleration of the universe. 




2 4 6 8 5 10 15 

T 

R 



Figure 1: Potential V(R) (left) and numerical solution to the scale factor R(t) (right) for a 
model with k — 1, £ = 0.1, A = 1 and p = 1. The five-dimensional Planck scale is taken 
k 5 = 1. 

Incidentally, when k 2 < 12 and «fAjj > 36£ 2 (A > 0) the universe starts at a big-bang 
singularity at r = and recollapses after reaching a maximum size of R max = [fi;fpo/(6£ + 
«||Ao|)]^. 

Let us analyze the case when k 2 > 12 which corresponds to a positive bulk cosmological 
constant. Now the asymptotic behavior of the potential is 

V{R) ~ I (3.18) 

and the situation is reversed with respect to the case analyzed before. For large values of R, 
the potential is always negative with V'(R) > and the decelerated evolution is controlled by 
the matter energy density through the p 2 term in the Friedmann equation (|3.4j) . If again we 
demand «|Aq < 36£ 2 (A < 0), we find that the universe does not have a big-bang singularity, 
but it starts at a finite size to undergo an inflationary expansion in which the scale factor 
grows as a power law. This regime is followed by a "graceful exit" phase where the universe 
enters a decelerated expansion which asymptotically leads to the p 2 term dominance. An 
example of this kind of solutions is presented in Fig. 2. In the left panel we have plotted the 
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potential V(R) and on the right one a numerical solution of Eq. ( |3.9| ) is shown. The total 
amount of inflation is just determined by the ratio between the size of the universe at the 
end of inflation V'(Rf) = and the original size of the universe V{Rj) = 0. 




1.5 2 2.5 3 3.5 4 5 10 15 

R T 



Figure 2: Profile of the potential V(R) (left) and numerical solution of the scale factor R(t) 
(right) for k = 4, £ = 1, Ao = — 1 and po = 3. We have also taken re 5 = 1. 

If instead we have k^Aq > 36£ 2 , the universe starts at R = and experiences a decelerated 
evolution to reach asymptotically the p 2 -dominated regime. This difference in the behavior 
when R — » can be understood by noticing that in the previous case the dynamics of the 
universe at small values of R is dominated by a large negative vacuum energy that prevents 
the formation of the singularity. 

Perhaps the most surprising feature of the solutions with k 2 > 12 is the fact that the 
late-time (large R) regime is dominated by the p 2 term in (|3.4j ). This might seem counter- 
intuitive, since it is usually argued that this term should only be important in the very 
early universe where the energy density of the universe is very large. The key element to 
understand the large- R behavior of these solutions is to realize that in our model the four- 
dimensional gravitational coupling evolves with time as R s , whereas the p term couples 
directly to the five- dimensional Planck scale, which is time independent. Thus if k is large 
enough, both the "quintessence" term and the term linear in the energy density become 
subleading with respect to the p 2 contribution when the scale factor grows. Therefore, 
the cosmological evolution is asymptotically dominated by the bulk term in the Friedmann 
equation. 
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Before finishing with the analysis of radiation branes, we study, for the sake of complete- 
ness, the marginal case k 2 = 12 where the bulk dilaton potential vanishes. As explained at 
the end of Section 2, when k 2 = 3m = 12 the brane tension is given by Eq. (|2.26|) , which 
contains a logarithmic dependence on R. Then V(R) is given by 



V(R) k *-_ 



12 



Z 2 -^^o + p -4 Po \ogRy 

3o 



R- 



(3.19) 



This implies that for both large and small R, the dynamics of the scale factor is controlled by 
the logarithmic term in the potential. The profile of the potential presents then two allowed 
regions for a zero-energy particle, for large and small R, separated by a barrier. 



4 Cosmologies on the brane in bulks with broken Z 2 - 
reflection symmetry 

As we discussed in Section 2, braneworlds can be also embedded in a bulk metric of the form 



( p. 2D without ^-reflection symmetry by considering different values of the energy scale £ at 
the two sides of the brane location. In this case, the evolution equations for the scale factor 
are obtained by solving Eq. Q2.17 ) for the Hubble parameter with the result 

= -^^R-^ k2 + |(A + pf + A (g _ e) 2 (A + P )- 2 R-P 2 . (4.1) 




Comparing with Eq. ( |3.1[ ), we notice the presence of an extra term which scales with a 
negative power of the brane energy density (cf. |r?fl). The expressions for p(R) and A(0) are 
the same as for the symmetric case. Proceeding as in the previous section, we reduce the 
problem of integrating the scale factor to finding the trajectory of a classical nonrelativistic 
particle with zero energy in the potential, 

a + e 



V(R) 



2 36 L 1 1 



Po 



R -(m-l) 



9 
44 



(£-£) 2 iH( fe2 - 3 ) [|A |iH (fc2 - 3) + Po R~ {m - 1) 



(4.2) 



where, as above, poR~ m represents the rescaled effective energy density on the brane. We 
notice that the last extra term in the potential associated with the breaking of reflection 
symmetry is always negative and therefore its contribution tends to increase the value of 
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. 2 

R 



H 2 = ( | ) . Moreover, it can be easily checked that this new term is bounded by the second 



one on the right-hand side of Eq. (|4.2| ) for all values of R provided |A | > Jr\/l£+ — £-1- 
As a consequence, when this bound is satisfied the dynamics of the brane universe is never 
dominated by the extra term in Eq. ( f4.1| ) and the cosmological evolution is qualitatively 
similar to the one found in the previous section. 

In order to show the qualitative differences with respect to the symmetric case studied in 
the previous section, we have to evaluate the relevance of the extra term in Eq. ( }4.1|) . Here 
we will focus on the solutions with phenomenological applications, therefore, for reasons 
that will be explained in Section 5, we restrict our attention to cases with k 2 < 3m. If 



|Ao| > Jta/|£+ — the extra term is negligible at small R compared with the dominating 
p 2 term. In general, however, accelerated expansion takes place for large values of R. In 
this latter case, the effective four-dimensional vacuum energy depends on the scale of Z2- 
symmetry breaking, (£? — £ 2 ) 2 , 



•> 36 A 4niXi 



V(R)~t L "™ J ^° lJ (4.3) 



■a 

In 

-4\2riQ/'e2 1 £-2 \ ,_4\2l ^ Oil' £2 £-2 \2 : 



Taking K5Aq[18(£ 2 + £ ) — k^o] < 81(£ 2 — £ 2 ) 2 it is easy to check that the universe rec- 
ollapses within finite proper time without passing through a phase of accelerated expansion. 
Therefore, these solutions are not very interesting from a phenomenological point of view. 
On the other hand, if k|Aq[18(£ 2 + £ 2 ) — n^X 2 ] > 81(£ 2 — £ 2 ) 2 , the universe starts again in a 
p 2 -dominated phase that, when p decreases, will be followed by a standard expansion domi- 
nated by the term linear in p in Eq. ( |4.1| ). Eventually, however, the (positive) quintessence 
energy density will dominate and the universe enters a regime of accelerated expansion. 
Actually, when R is large enough and higher powers of p can be neglected [which in- 

3 

eludes the limit R 3> (po/|Ao|) Sm ~ k2 } the dynamics of the universe can be derived from the 
approximate Friedmann equation, 

2 

« ^-G N (R)p Q (R) + ^-G N (R)p(R), (4.4) 




where the new Newton constant is given by 

^5 | Aq | 



Gn(R) 



48tt 



(£i-£ 



2 \2 



R~^ (4.5) 
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and the quintessence component of the energy density is 



pq(R) 



2«||Ao| 



K 5 A V K 5 A / 



R-~. (4.6) 



For consistency, only terms linear in ^p^-(£ + — £_) 2 are retained upon substitution in Eq. 
( |4~4"D . We see that both Newton's constant and the quintessence energy density reduce to 
expressions Q3.5|) and ( |3.7| ) for £ + = 

Let us consider now the case in which the difference between the cosmological constant 
at the two sides of the brane is "large", i.e. «;||Ao| < 3a/|£+ — Now the last term in 
Eq. (|4. 1|) will dominate for a certain range of values of R. Thus, the p 2 -dominated epoch at 
small R will be followed by a regime where the expansion is driven by the term proportional 
to (A + p)~ 2 . This will happen for the following range of values of the scale factor: 



2 \ 3m-fc 2 

^ ] <R< t^t) ■ (4-7) 



Ar 



For larger values of R, the universe will enter a quintessence-dominated regime. 

To summarize, we have seen that when the difference between the values of £ at both 
sides of the brane is small (in the sense explained above), the cosmological dynamics of the 
universe does not show a radically different behavior with respect to the symmetric models. 
Actually, for large values of R the effect of the breaking of the reflection symmetry shows in 
a correction of both Newton's constant and pq, although their scalings with the scale factor 
are not affected. On the other hand, for larger values of |£+ — the main difference is the 
existence of an intermediate phase interpolating between the p 2 regime and the "standard" 
expansion dominated by the linear term in the Friedmann equation, where the dynamics of 



the universe is controlled by the last term in Eq. (fO). At large values of R, however, the 



dynamics can still be described by the approximate Friedmann equation ( |4.4|) . 

5 Phenomenology of quintessential brane cosmologies 

After having investigated their mathematical properties, in the following we will study phe- 
nomenological aspects of the quintessential brane cosmological models presented in Section 3, 
where the acceleration of the universe at late times is induced by the rolling of the braneworld 
in a static bulk geometry. We will focus our attention on the ^-reflection symmetric models. 
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We saw that whenever k 4 Aq > 36£ 2 and k 2 < 3m (with m > 0), solutions describ- 
ing a flat Friedmann-Robertson- Walker universe can be found which, after going through 
radiation/matter-dominated epochs, undergo accelerated expansion at late times. For these 
types of solutions, the energy content of the very early universe is dominated by radiation 
and the expansion is controlled by the p 2 term in Eq. (|3.4j) . When the density is low enough, 
a phase of "standard" expansion starts where the dynamics of the universe is governed by 
the term linear in p in the Friedmann equation. It is during this phase that nucleosynthesis 
has to take place, as well as the decoupling of matter from radiation. Finally, some time after 
the universe enters in a matter dominated epoch, the quintessence term starts dominating 
the dynamics and the universe begins its accelerated expansion. This acceleration is driven 
by the bulk scalar field through the term pq in the Friedmann equation, which, as we will see 
below, plays the role of effective quintessence on the brane. In addition to this, the effective 



Newton constant in our models, Eq. ( |3.7D , evolves with time through its dependence on the 
scale factor. 

We normalize the scale factor on the brane so it takes the value Rq = 1 at the present time 
and we use the temperature of the radiation to label the different epochs in the evolution of 
the universe, 

T(R) = !™ (5.1) 

From our results in Sec. 2.3 we saw that the energy density on the brane scales with the 
temperature as 



T 



3(l+«j e ff) 



where w e g is give by Eq. ( |2.24| ). As discussed above, for the particular case of radiation w c s = 
w — |, and we recover Boltzmann's law, p{T) ~ T 4 . However, when w = the effective 
barotropic index is nonvanishing. Nonetheless, in the following we will denominate "matter" 
to any perfect fluid with w e s = 0, whose energy therefore scales with the temperature as 
p(T) ~ T 3 , as seen by a brane observer. Assuming by naturality that b ~ k, we find that 
w e g — wc^k 2 . As we will see later, this number is of order < 1CT 2 so the difference between 
the "microscopic" barotropic index and the effective one (which governs the cosmological 
evolution) is very small. 

To study the evolution of the universe in more quantitative terms, it is useful to introduce 
three parameters, 6q(T), Q p (T), and 6 p 2(T), which measure the contribution to the expan- 
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sion of the three terms on the right hand side of the Friedmann equation ( |3.4| ) , normalized 
to the Hubble parameter, 




G N (T)p Q (T), 



G N (T)p(T), 



(5.3) 



6 p2 (T) 



36H 2 



where Pq(T) is given by Eq.( |3.5| ). In this way the Friedmann equation is written as 



The crossover between the p 2 -dominated epoch and the standard regime will take place 
at a temperature T c at which 6 p 2(T c ) ~ p (T c ). By using the definitions ( |5.3|) we find 



On the other hand the standard expansion of the universe will be terminated at a temperature 
Tq < T c in which 6 p (Tq) ~ Qq(Tq). Since this has to occur after the universe has entered 
a matter-dominated era, we find 



As a consequence, the period of standard expansion can be very large if the net cosmological 
constant on the brane, k 4 Aq — 36£ 2 , is small enough. 

In order to illustrate the discussion above, in Fig. 3 the evolution of the Hubble parameter 
is shown for some "realistic" values of the parameters in a cold dark matter model. We can 
see clearly several regimes. At high temperatures (T > 10 12 K), the expansion is dominated 
by the p 2 term in the Friedmann equation. This is followed by a long period of standard 
expansion for 10 12 K < T < 10 K. In this regime, radiation (m = 4) dominates the expansion 
for temperatures higher than T eq ~ 10 4 K, whereas for lower temperatures the energy density 
of the universe is dominated by matter (m = 3). Finally, when T < 10 K, the quintessence 
component takes over and the universe starts its accelerated expansion. 



i = e Q (T) + e p (T) + e p2 (r). 



(5.4) 




(5.5) 




(5.6) 
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Figure 3: Evolution of the Hubble parameter (GeV) vs. temperature (K) for a realistic 
quintessence model. The values of the parameters are: p — 0.43pg )0 = 2.43 x 1CT 47 GeV , 
G Nj0 = 6.71 x l(T 39 GeV~ 2 , T eq = 10 4 K, k 2 = 0.01 and k 5 = 10- 8 GeV"t. 

5.1 Bounds from big-bang nucleosynthesis 

The abundances with which light elements are synthesized in the early universe are extremely 
sensitive to changes in the expansion rate of the universe and/or in the fundamental constants 



2% piql . Thus, primordial nucleosynthesis provides a good way to bound the values of the 
constants and scales that characterize a given cosmological model. In the case at hand, the 
effective gravitational coupling varies with time and from Eq. ( |3.7|) we can see that this 
change is determined solely by k, namely 

p- = ~H. (5.7) 

In order to explain the primordial abundances of 4 He and other light elements, however, 
Newton's constant at the time of nucleosynthesis cannot be very different from the present 
value. Actually, its time variation has to satisfy [[2(J 




< ±0.01 h~ L H . (5. 
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Therefore, we find that the parameter k that determines the slope of the potential of the 
bulk dilaton has to be rather small, 



k 2 < 0.03 h" 1 , (5.9) 



which, from Eq. ( |3.6|) , implies Wq + 1 < 0.003 h" 1 (cf. 

A second condition to be met by a cosmological model in order to account for the pri- 
mordial abundances of the light elements is that the energy density of the universe at the 
time of nucleosynthesis has to be dominated by radiation. In particular, the dark energy 



contribution to the density parameter satisfies the following bound |28[ (see also |29[ and 
references therein): 

^(TBBN)dark < 0.045 (at 2<j confidence level), (5.10) 

where Tbbn — 10 10 K. In our model, the dark energy component corresponds to the contri- 
bution of pq. Thus, we define 

^=^m' (6 - n) 

where p{T) is the matter/radiation energy density. 

From Eq. (|5.11|) , and assuming a cold dark matter scenario, we can obtain an expression 
of ^q(Tbbn) in terms of quantities at the present time as 

3 4 ^fc 2 " 1 " 1 

1 — Qq(T ) /T cq \ /Tbbn \ /Tbbn 



^q(Ti 



BBN 



1 



^q(To) V To / V T cq / \ T c 



(5.12) 



where T eq ~ 10 4 K is the temperature at which the density of radiation and matter are equal, 
and T = 2.73 K is the temperature of the CMB at the present time. Since the measurement 
of CMB anisotropies and the supernova data indicate that Qq(Tq) ~ 0.7, we find that for 



models satisfying the bound Q5.9D the contribution of the dark energy component at the time 
of nucleosynthesis is 

^q(Tbbn) ^ 10" 35 , (5.13) 
which falls well inside the nucleosynthesis bound ( [5.10| ). 
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5.2 Bounds from observations 

We have seen that by satisfying the bound ( |5.9|) , our quintessential brane cosmologies com- 
ply with the constraints imposed by big-bang nucleosynthesis. It is interesting to notice, 
however, that none of the two bounds considered in the previous section involves any other 
parameter of the model but the numerical constant k, or equivalently the quintessence effec- 
tive barotropic index wq. In particular we do not get any information as to the size of the 
dimensionful parameters Ao and £ or the five-dimensional Planck scale k§. 

In order to extract bounds for these scales, we have to compare the values of pq and 
the present value of the four- dimensional Newton constant with the observed values. From 
Eq. (|3.7|) we get the value of the Newton's constant by setting R = 1. Using Gn,o = 



6.71 x l(T 39 GeV~ 2 , we find 

4\X \ = Qk 2 4 ~ 1(T 36 GeV" 2 . (5.14) 

The second observational result we can use to estimate the scales of the model is the 
value of the dark energy component of the density parameter, fl(T ) dark ~ 0.7. From Eqs. 
( CT ) and (El), we find that 



~4.67«t|Ao|po, ( 5 - 15 ) 

where po is the matter density of the universe at the present time. Using Eq. (|5.14|) and the 
value po ^ 2.43 x l(T 47 /i 2 GeV 4 |1 we arrive at 



K 4 A 2 -36£ 2 ~l(r 82 GeV 2 . (5.16) 

Incidentally, using Eqs. Q5.15| ) and ( |5.6| ) we find that Tq ~ 3.63 K. 

Equations ( |5.14|) and Q5.16| ) are the only two estimates that we can get just using "in- 
brane" information. In order to go forward and get some estimate for the five- dimensional 
Planck scale K5 or the brane tension Ao, we need some knowledge of the fundamental energy 
scale in the higher-dimensional theory. The p 2 term in the Friedmann equation (|3.4| ) provides 
a window to the bulk dynamics, since it couples directly to the five- dimensional Planck scale 
K5. Therefore, it would be very interesting to have some information about the transition 
between the p 2 -dominated epoch and the standard expansion. Even if we do not have a 
concrete value for the temperature at which this transition takes places, we can still get a 
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bound for Ao by going back to Eq. (5.5|). In order to preserve big-bang nucleosynthesis, 
the p 2 -component in the Friedmann equation has to be negligible at Tbbn p0 |- Imposing 
T c > 10 10 K, we are led to 

|A | > 10 40 p ^ KT 7 GeV 4 . (5.17) 



Using Eq. (|5.14 ) the following bound for the five- dimensional Planck mass is found, 

M 5 = k~* > 10 5 GeV, (5.18) 
and consequently £ ^> 1CT 22 GeV. Notice that this bound for the five-dimensional Planck 



mass is stronger than the one obtained in the quintessence scenario of Ref. |!l|, although it 
is again weaker than the one obtained from the Randall-Sundrum model ||. 

6 Conclusions 



In the present paper, we have studied the bulk solutions constructed in Ref. JL4] in a 
different physical setup. Instead of being at a fixed location, the braneworlds move in a five- 
dimensional dilatonic bulk space time with and without Z2-reflection symmetry. As a result, 



and contrary to the case analyzed in Ref. |TJL the cosmological dynamics is not driven 



by a time-dependent background but by the movement of the braneworld in an otherwise 
static bulk space-time. The resulting class of brane cosmologies admits, depending on the 
choice of parameters, either inflationary behavior at early times and decelerated expansion 
at late times or a decelerated expansion at early times which evolves into a quintessencelike 
behavior at late, i.e., present, times. 

The realization of quintessence in these models seems to be particularly interesting since it 
does not rely on the introduction of an additional scalar field on the brane. It is in fact driven 
by the bulk dilaton. As the brane moves along in the static bulk, it sweeps out sections with 
different values for the bulk dilaton and, from the point of view of an observer living on the 
brane, this is interpreted as a time- dependent cosmological constant. For a suitable choice of 
parameters, this evolution is such that the energy density of this effective quintessence only 
becomes important at the present time. However, in addition to the induced time- dependent 
cosmological constant the bulk dilaton also causes a time-dependent gravitational coupling. 
Using observational constraints on the variation of Newton's constant and the present value 
of the cosmological constant gives bounds on the parameter of the brane cosmologies at hand. 
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The fact that k|Aq — 36£ 2 <C k^Xq,^, 2 seems to indicate that the nonvanishing value of the 
four- dimensional vacuum energy is the result of a slightly broken symmetry Furthermore, 
we find a lower limit for the five-dimensional Planck mass of 10 5 GeV. 

Here we are taking a phenomenological approach and therefore we do not attempt to 
explain these values for the scales of the model from first principles. However, it would be 
very interesting to investigate whether the hierarchy of scales found here can be obtained in 
a natural way by embedding these five-dimensional models in a string/M-theory setup. 

We have shown that for a certain range of parameters, inflation takes place in the early 
universe. However, this phase is followed, after a graceful exit period, by a decelerated 
Friedmann- Robertson- Walker universe in which the expansion is driven by p 2 . This fact 
makes this particular set of solutions not very promising as phenomenological models of the 

0- 



universe 



We have also considered brane cosmologies with broken Z 2 -reflection symmetry and found 
that for a slight symmetry breaking the evolution follows qualitatively that of the symmetric 
case. In this case, for large values of the scale factor the evolution is approximately described 
by a Friedmann equation where both Newton's constant and the quintessence energy density 
receive corrections depending on the scale of Z 2 -symmetry breaking. 



One of the interesting properties of the braneworlds investigated in [B3 is the existence 



of an extension of the mechanism of self-tuning [19j in which the matching conditions at the 
brane location force the cancellation of all vacuum energy terms on the right-hand side of the 
Einstein equations in four dimensions. However, from our analysis in Section 3 we find that 
the physical mechanism enforcing this self-tuning is the condition that the braneworld is at 
a fixed position in the bulk. Indeed, from Eq. ( |3.1|) we conclude that the condition R = 
for any R is equivalent to consider a vacuum brane (w = —1) with zero net cosmo logical 
constant in four dimensions, k^Xq = 36£ 2 . From this point of view, at least in the case of a 
static bulk geometry, the self-tuning can be seen as resulting from the tuning of the velocity 
of the brane in the bulk to zeroQ However, this is not necessarily a fine-tuning, since the 
condition R = can be automatically implemented if the braneworld is confined to lie on 
the fixed points of an S 1 /^ orbifold, as it is the case in Hofava-Witten type models |13 . 



5 This result extends to the cases with Z2-reflection symmetry breaking, where again a nonmoving brane 
is only achieved by considering a vacuum brane with k|Aq = 9(£+ ± £-) 2 - 
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